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The one-dimensional spin-S Ji — J2 XY model is studied within the bosoniza- 
tion approach. Around the two limits (J2/J1 ^ l,J2/Ji ^ 1) where a field 
theoretical analysis can be derived, we discuss the phases as well as the dif- 
ferent phase transitions that occur in the model. In particular, it is found 
that the chiral critical spin nematic phase, first discovered by Nersesyan et al. 
(Phys. Rev. Lett. 81, 910 (1998)) for S=l/2, exists in the general spin-S 
case. The nature of the effective field theory that describes the transition be- 
tween this chiral critical phase and a chiral gapped phase is also determined. 



PACS No: 75.10.Jm 



I. INTRODUCTION 



The interplay between frustration and quantum fluctuations in low-dimensional spin systems 
has attracted much interest. One of the main reasons for this attention is that frustration 
is expected to lead to new exotic phases as well as unconventional spin excitations. In the 
one dimensional case, powerful non-perturbative methods are available and the key features 
of frustration can then be analysed in depth. In this respect, the phase diagram of the one 
dimensional S=l/2 Ji — J2 XY model has been studied extensively over the years. The lattice 
Hamiltonian of this model is defined by: 

n n 

where ^ S^± iSl is a spin-1/2 operator at site n and J2 is a competitive antiferromagnetic 
interaction (Ji, J2 > 0) which introduces frustration in the model. The phase diagram of 
the Hamiltonian (|^) is expected to be rich. For a small value of J2, one has a spin fluid 
phase (XY phase) characterized by gapless excitations with central charge c = 1 whereas for 
J2/ Ji ~ 0.32 a phase transition of Kosterlitz-Thouless (KT) type occurs and the model enters 
a massive region with a two- fold degenerate ground state (dimerized phase) [1). Interestingly 
enough, it was recently predicted by Nersesyan et al. within the bosonization approach 
that in the large J2 limit, where the model can be viewed as a two-leg XY zigzag ladder, a 
critical spin nematic phase with chiral long-range order A Sn+i)z) 7^ 0) should emerge. 

This unconventional phase with unbroken time reversal symmetry is characterized by nonzero 
local spin currents polarized along the anisotropy z-axis. The transverse spin-spin correlation 
functions are incommensurate and fall off with the distance as a power law with the exponent 
1/4. However, this phase has not been reported in the numerical calculations of Rcfs. |^,^. 
In contrast, only the previous well known phases (the spin fluid and dimerized phases) have 
been found. A numerical analysis of the model using an exact diagonalization method with 
twisted boundary conditions has pointed out the presence of incommensuration in the large 
J2 limit [0 but within this approach it has not been possible to conclude on the criticality or 
not of this phase. Very recently, Nishiyama |^ has investigated the existence of chiral order of 
the Josephson-junction ladder with half a flux quantum per plaquette by means of the exact 
diagonalization method. He was able to show that the critical phase predicted in Ref. |^ does 
exist in the range of the parameters of the model (|^) in constrast to the numerical findings of 
Ref. |. 

The situation is less controversial in the S=l case and the corresponding phase diagram has 
been determined by a DMRG study The model with J2 = is a critical spin fluid (the 

so-called XYl phase |^) and as soon as the next-nearest-neighbor interaction is switched on 
the Haldane phase is stabilized (KT transition). As noted in Ref. ||], this fact seems to be in 
contradiction with the bosonization result obtained in the small J2 limit which suggests that 
the XYl phase extends to finite J2 / Ji . Increasing on the value of J2/ Ji , the authors of Refs. 
1^,^ have reported the occurence of two successive transitions: A first one at {J2/ Jijci — 0.473 
(Ising transition) from the Haldane phase to a gapped phase with chiral long-range order (chiral 
gapped phase) and a second transition at {J2/ Ji)c2 — 0-49 (presumably a KT transition) with 
a chiral critical phase which corresponds to the spin nematic phase discussed by Nersesyan et 
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al. [| in the context of the two-leg S=l/2 XY zigzag ladder. In this latter phase, the spin-spin 
correlation are incommensurate and decay with a power law with an exponent approximatively 
equals to 0.15 ||]. Recently, the phase diagram of the model (|^) in the general spin-S case has 
been further discussed in Ref. by means of a large-S approach. It has been found that the 
existence of the gapless and gapped chiral phases is not specific of S—1 but is rather a generic 
large-S feature. The predicted phase diagram has four different phases. First of all, one has 
the XYl critical spin fluid phase and for a finite value of J2 (KT transition) one enters the 
Haldane phase. As J2/J1 is further increased, there are the two successive phase transitions 
previously discussed: An Ising-type transition separating the Haldane phase from the chiral 
gapped phase, followed by a KT transition corresponding to the transition between the chiral 
gapless phase and the chiral gapped phase. However, this large-S approach does not distinguish 
between integer and half-integer spins and also it does not take into account the possibility 
of a spontaneously dimerized phase. In addition, for the very special S=l case, the numerical 
calculations of Refs. predict rather that the critical phase XYl does not extend for a finite 
value of J2. 

In this paper, we shall investigate the phase diagram of the one-dimensional spin-S Ji — J2 XY 
model within the bosonization approach. Using the Abelian bosonization of a general spin-S 
operator introduced by Schulz , the low-energy physics of the Hamiltonian (|l|) can be studied 
in two different limits: in the weak coupling limit when J2/J1 ^ 1 and in the ladder limit 
J2/J1 ^ 1- This enables us to determine the nature of the phase transitions that occur in the 
model from a field theoretical point of view together with a comparison with the numerical 
results as well as the large-S preditions |l^. The remainder of the paper is organized 
as follows: the weak coupling analysis (J2/J1 ^ 1) is given in Section II whereas the zigzag 
limit of the model (J2/J1 ^ 1) is performed in Section III. Section IV presents our concluding 
remarks and finally the conventions and some technical details used in this work are described 
in the two appendices. 



II. WEAK COUPLING LIMIT 



In this section, we shall investigate the low energy physics of the one-dimensional spin-S 
Ji — J2 XY model in the limit J2/J1 ^ 1 within the bosonization approach. This enables us 
to study the stability of the critical XYl phase (for J2 = 0) upon switching on a small next- 
nearest- neighbor interaction. The Hamiltonian (|^) in the S=l case can be tackled within the 
bosonization framework by representing spin-1 operators as a sum of two spin-1/2 operators: 
Si^ = s^^ + S2n- It h^-s been argued Jl2|-p^ that the additional local singlets introduced will 
lead to extra levels with higher energy than the triplet states so that the ground state and the 
low-lying excitations are correctly captured. In particular, with this representation, Timonen 
and Luther and Schulz predicted the correct phase diagram of the one-dimensional 
anisotropic antiferromagnetic spin-1 Heisenberg model with single-ion anisotropy. Moreover, the 
effect of weak randomness on this latter model has also been analysed within this bosonization 
approach p^ . This procedure was further generalized by Schulz by representing a general 
spin-S operator as the sum of 2S spin-1/2 operators s^„, a — 1, .., 25*: 
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25 
a=l 

from which the low-energy physics of the one-dimensional spin-S Heisenberg model can be 
captured and in particular the difference between half- integer and integer spins [p7| . 

Let us first review the continuum description of the spin-S XY chain first obtained by Schulz 
Using the decomposition (||), the Hamiltonian for J2 = writes: 

25 

nXY, = Y E E + H.C.) + f E E {<n\n+l + H.c) . (3) 

a— 1 n a^b n 

The first term in this equation corresponds to 25" decoupled spin-1/2 XY chains. As recalled 
in the Appendix A, the low-energy physics of the spin-1/2 XY chain can be extracted from the 
introduction of a single U(l) bosonic field ip with chiral components fR.L- As a consequence, the 
next step of the approach is to introduce 2S chiral decoupled bosonic fields ifaR.L, a = 1, •■7 25' 
and the Hamiltonian density of the first term in Eq. (^) in the continuum limit is given by: 

2S 

Hxy^yT. ((^-'/'«)' + (d.^af) , (4) 

a=l 

where = ^aL + VaR., '&a = 'PaL - VaR bciug the dual field, and Vo = JiOo (oq being the 
lattice spacing) is the spin velocity. One can then derive the continuum limit of the second 
term in Eq. using the bosonic description (A7) of the spin-1/2 operators described in 
the Appendix A: 

+ —^^= (exp (iSy^ifaL + i\/7r(^a_R) + exp {~i3y/TTipaR - i\/7r.^aL)) ■ (5) 

Notice that a priori this procedure has only a sense provided that the coupling constant asso- 
ciated to the second piece of Eq. (||) is much smaller than Ji and this is clearly not the case 
here. However, one expects in this problem, on general grounds, a continuity between weak 
and strong coupling limits so that it is natural to bosonize the second term of Eq. (||) using 
the correspondence (H). The leading part of the density Hamiltonian associated to (||) in the 
continuum limit reads thus as follows: 

2S , 

UXYI - y E ((^-^a)' + {^-^af) - ^ E '^O^ " ^^^) ' 

a=l ° ai^h 

It is then suitable to switch to a basis to single out the degrees of freedom that will remain 
critical in the infrared limit. To this end, let us introduce a diagonal bosonic field and 
25—1 relative bosonic fields '^mR(L)^ '^i — (1, 25 — 1) as follows: 

^^R{L) = ^ (^1 + - + ¥'2s)ii(i) 

1 

The transformation is canonical and preserves the bosonic commutation relations. This 
basis has been introduced in Ref. [|l8| in the Abelian bosonization study of the one-dimensional 
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Hubbard model with a SU(N) symmetry. The inverse transformation of Eq. is easily found 
to be: 



1 $ 



W.(L) 



1 ^ 

^aR{L) = —J^^+RiL) - \ *(a-l)_R(L) 



2S 

2S— 1 

I 



1 /2^-l^ /o^ 

V2Sfl(L) - ^ Y 25' '^(2S-l)R.(L)- (8) 

Using Eq. one observes that all relative dual fields Om, m — I, .., 25— 1 arc pinned whereas 
the diagonal bosonic field (<i>+ = + ^+r) is a strongly fluctuating field so that: 

HxYi ^ ^ {{d.^+f + (d.e+f) , (9) 

where 0+ = <i>+_L — ^+r- From Eq. (||) and by integrating out the massive degrees of freedom, 
the expression of the effective spin-S density in terms of the massless bosonic field in the + 
sector can be deduced: 



(-f exp (±iv/V25 e+) . (10) 



The dual field 0+ is thus a compactificd bosonic field with radius Rs = \j2SpK. Using the 
general relation between the radius (2t:RsRs — 1), we deduce that the compactified radius 



of the bosonic field $+ is: Rs = I/VS-kS. Futhermore, one deduces from Eq. (10) that the 
transverse spin-spin correlation function has a power law behavior with an exponent r]± = 
1/(45). The value of this exponent coincides with the prediction of Alcaraz and Moreo (Toj 
who have analysed the critical properties of the XXZ spin-S Heisenberg model by means of a 
combination of conformal invariance and exact diagonalizations techniques. It is worth noting 
that the value of the exponent in the S~l case {ri± — 1/4) has been predicted by Kitazawa 
et al. within a level spectroscopy analysis of the S=l bond-alternating XXZ spin chain. 
Finally, one should observe that the uniform part of the spin density is a short-ranged piece 
since the massive modes that enter in this expression have a zero vacuum expectation value and 
thus give rise to an exponential decay in the uniform part of the spin-spin correlation function. 
However, as shown by Schulz in the special half-integer case, higher orders of perturbation 
theory produce a strongly fluctuating piece in the uniform part of the spin density ( pO| ) with 
scaling dimension 25 -I- 1/(85) which is less relevant than the alternating contribution in Eq. 
( p^ which has scaling dimension 1/(85). 

With all these results at hands, the stability of this critical XYl phase with respect to a 
next-nearest-neighbor exchange interaction J2/J1 ^ 1 can be analysed. To this end, let us first 
rewrite the second term {Hxy2) of the Hamiltonian (|l|) in terms of the 2S spin-1/2 operators: 

2S 



nxY2 - Y E E i<nS-n+2 + H.c) + {<n\n+2 + H-c) . (11) 
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The first part of tliis equation corresponds to tlie sum of N decoupled next-nearest-neiglibor 
8=1/2 XY chains. The resulting continuum limit has been obtained by Haldane in the erratum 
of Ref. Q and is reviewed for completeness in the Appendix B (see in particular Eq. (B4)). The 
continuum limit of the second term in Eq. (11) can be obtained using the bosonized description 
(^) of the spin density sj. The resulting continuum limit of the Hamiltonian ( pT| ) reads thus 
as follows: 

HXY2 ^ V [d^daf + A V cos [V^ {da - M) 



nan 

a=l ^' a<b 



J 25 
•J2 \ ^ 



^2 / . '^Q^ ISTTy^g) + cos ( VItt {ipa + ^pb) \ COS ( Vtt (^9a - Z^j,)) . (12) 

Using the canonical transformation (||), one finally obtains the following effective Hamiltonian: 

H^^(^i^(9.e+)V-i(a.$+)'), (13) 

with 



8Js 
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Of course, as usual, these latter identifications hold only in the vinicity of the gaussian fixed 
point at J2I J\ = 0. However, one should carefully look at the higher order corrections of 
perturbation theory that might generate an additional operator in the effective theory (|l^) and 
potentially destablizes the XYl phase. In this respect, integer and half-integer spins should 
be treated separately. For integer spins S, the last term in Eq. ( p^ ) leads to the following 
contribution at the Sth order of perturbation theory: 

» s s 

/ Y\ d^^i W cos {(P2i-1 + </52i)) {Xi) COS ( (t92i-l - d2i)) [xi) , (15) 

•' i=l i=l 

which, after integrating out the short-ranged degrees of freedom and using the canonical trans- 
formation (^), gives rise to a fluctuating field in the -|- channel: cos (v^8^$+). In the same way, 
for half-integer spins S, one has the following contribution at the 2Sth order of perturbation 
theory: 

. 2S 2S 

/ JJ^d^Xi J|cos (yiTr{ipi + ipi+i)^ {xi)cos (i9i - t?i+i)) (xi) , (16) 
i=i i=i 

with the identification: ip2S+i = and 'd2S+i = '^i- One then obtains the following operator: 
cos (\/32^$+) after averaging on the short-ranged degrees of freedom. The effective field 
theory associated to the spin-S Ji — J2 XY chain in the weak coupling limit is thus: 

n^^-(^K [d^Q+f + {d^.^+f^ - ^ cos (/iV8^$+) , (17) 

with fJ, — 1 (respectively 2) if S is integer (respectively half-integer). One should note that 
the Hamiltonian (|l7| ) corresponds to the effective field theory of the spin-S XXZ Heisenberg 
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chain derived by Schulz f^. In fact, the last operator in Eq. ( [l7|) can also be justified from 
a symmetry analysis. Indeed, under the one-step translation symmetry, the bosonic field ipa 
transforms according to (see Eq. ( |A8D of the Appendix A): (fa ^ y^a + v^/^ +PaV^j Pa being 
integer. From the definition (|7|) of the diagonal bosonic field <I>_|-, one thus has: 

1>.-'J>. + /f + (18) 

from which we conclude that the cos ^/i\/87r5$+^ term is the operator invariant under the 
translation symmetry with the smallest scaling dimension. 

The phase diagram of the spin-S Ji — J2 XY chain in the small J2/J1 limit can then be 
deduced from the structure of the effective field theory (|l7|). For a small value of J2/J1, the 
cosine operator in Eq. (|l^) is a strongly irrelevant contribution and the system is critical with 
central charge c = 1 (Luttinger liquid): it is the spin fluid XYl phase that extends to a flnite 
value of J2. As Jij J\ increases, one expects from Eqs. ( p^ , pT[ ) a KT phase transition from this 
spin fluid phase to a fully massive region (dimerized or Haldane phases depending on the nature 
of the spin S [^). At the transition, the Luttinger parameter is equal to: = l/{Sfi^) 
and a very naive estimate of the critical value of {J2/ Ji)c can then be deduced from Eq. ( p^ ) 
within the bosonization approach: 

A^.'C^V-l)^ (19, 



In particular, for 8=1/2, one finds {J2/Ji)c = 37r/32 ~ 0.2945 which is not too bad in com- 
paraison to the value obtained in the numerical simulations of Ref. ||l|: {J2/Ji)c — 0.3238. 
Moreover, in the S=l case, the XYl phase is destabilized upon switching on a nonzero value 
of J2 in full agreement with the numerical findings of Refs. 1^,^. The origin of the discrepancy 
noted in Ref. ||^ between the DMRG study and the bosonization results obtained in Ref. 
[]lO|| stems from the fact that the latter authors do not look at higher orders in perturbation 
theory as in this work. The S=l case does not correspond to the generic situation since we 
observe from the estimate ( pj|) that the size of the XYl phase increases as S increases in the 
half-integer and integer cases. In this respect, one should note that the situation is in close 
parallel to the phase transition between the XYl and the Haldane phases in the integer spin-S 
XXZ Heisenberg chain. In the S=l case, the resulting phase transition occurs precisely at the 
XYl point P, ^ , pO| whereas the XYl phase extends considerably as S increases [p|, p^p^ , |2^ . 



III. THE ZIGZAG LADDER LIMIT 



We shall now study the model (gj) in the ladder limit Ji <C J2 where it can be viewed as a 
two-leg spin-S XY ladder coupled in a zigzag way. For S=l/2 Heisenberg spins, the effect of a 
transverse zigzag interchain interaction has been extensively studied in Refs. [^ , p6| j3| j2^ , ^^ and 
also in Ref. |^ in the S=l case. In the special case of S=l/2 XY spins, it has been found by 
Nersesyan et al. [|| that the model is a critical spin nematic. In this section, we shall investigate 
the existence of such a phase in the general spin-S case and study its stability as the interchain 
interaction is further varied. 
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A. critical spin nematic phase 



The lattice Hamiltonian of the model (|l|), considered as a two-leg spin ladder, is defined now 
as follows: 

J2 



^ - Y X! ('^l,n'^l,n+l + '^2,n-l/2'^2,n+l/2 + ^ 

E {^in + ^2:„+l/2) + H .c) , (20) 



Jl 

' 2 



where S^^ (respectively S^^_^^^r^ is the spin-S operator of chain of index 1 (respectively 2) at 
site n (respectively n + 1/2). It is more suitable to change the labeling of the second chain in 
the following way to perform the continuum limit of the model: 

'^=iY.Y. {SlnS-,n+l + H.C) +^Y. (^U i^ln + ^2:„-l) + H .c) . (21) 
a— 1 n n 

At this point, one should note that the interchain zigzag coupling can also be written as (using 
intrachain periodic boundary conditions): 

-H^nt = Y E ( {^in + Sln+l) S^.n + H .c) . (22) 
n 

Consequently, we shall thus write the interacting part of the Hamiltonian ( pl|) in a symmetrized 
way for taking the continuum limit of the model: 

^= YEE(^ln^^n+l+^-^-) 

a— 1 n 

+ Y ^ ( (^^." + ^^'"+0 + ^1'" ^^2:„ + ^2:„-i) + H.c) . (23) 

n 

In the absence of the interchain coupling ( Ji = 0), the model corresponds to two decoupled spin- 
S XY chains. As seen in section II, it is critical with central charge c = 2 and its low-energy 
physics can be obtained with the introduction of two decoupled chiral gapless bosonic fields 
^a+R,L (a=l,2). The leading contribution of the spin density comes from the alternating 
part (see Eq. (|l^): 

St ^ i-lf'''" exp (±iV^ Qa+) , (24) 

A being a non-universal constant. From Eq. (p^, we deduce the continuum limit of the model 
( ^ ) in the small Ji ^ J 2 limit: 

^ - i E ((^-*-)' + (5-0a)') + 9 a.e+ sin (^1©-) ' (25) 

where g = JiA^-^/tt/ (45*) and we have introduced the symmetric and antisymmetric combina- 
tions of the two bosonic fields: 

^± = -^($l+±$2+) 

e± = (61+ ± 82+) . (26) 
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The Hamiltonian ( |25| ) describes a nontrivial field theory since the field with coupling constant 
g, called twist term in Ref. is a parity symmetry breaking perturbation with a nonzero 
conformal spin (equal to one) . The effect of such term is rather unclear since the usual irrelevant 
versus relevant criterion does not hold for such a nonscalar perturbation (see for instance Ref. 
Q). The simplest spin-1 conformal perturbation is the uniform part of the spin density {d^^) 
that couples to a uniform magnetic field along the z-axis. In this case, this term leads to 
incommensuration as is well known. It is thus natural to expect some incommensurability effect 
in the model (|5|) due to the twist term as emphasized by Nersesyan et al. In particular, 
the presence of incommensuration in the system can be found by a direct mean-field analysis 
of the model (p5|). Indeed, it is easy to see that the mean- field Hamiltonian separates into two 
commuting parts: TLmf = + with 



n. = \ -f (a.e_)^j - ^ sin ( ^^e_ ) , (27) 



the mean-field parameters being: 



.g(sin ( 



^ = -5(9.e+). (28) 

The Hamiltonian (7i+) is easily solved by the redefinition G)+ 8+ — kx/v. The -I- sector 
displays thus criticality with a nonzero topological spin current in the ground state: {dxQ+) = 
— k/v ^ 0. In contrast, the Hamiltonian (7i-) in the other sector is a standard sine Gordon 
model at f3'^ = tt/S which describes a massive theory with massive quantum solitons and 
their bound states (breathers) together with massive kinks. The dual field 9_ is locked at: 
(e_) = V7iW4 sgn^ (mod \/4SV). The mean-field analysis can then be closed using the fact 
that: (sin(A/7r/5' 6_)) — c(ao|/i|/w)^/'^'^~^-' (c being a constant that can be determined Q) 
and one easily finds: 

ao \ V / 

K = ±— — c5^. (29) 

ao V w / 

From the correspondence (|2^) , one can then estimate the asymptotic behavior of the transverse 
spin-spin correlation functions of the model which display an incommensurate critical behavior: 

„iqsx 

{Slix)S-m^j-^^j^, a =1,2, (30) 

with qs ~ Tr/ao ^ ( Ji/J2)^'^/*-''"^^^-'. The transverse spin-spin correlation functions fall off thus 
with the distance as a power law with the exponent 1/(85'). In the S=l case, one should note 
that this exponent (1/8 = 0.125) found in this bosonization study is in good agreement with 
the numerical findings 0.15 of the DMRG analyis of Ref. 1^. 

Besides this incommensurate critical behavior observed in the spin-spin correlation functions 
( |30| ) , the physical picture of this phase obtained at the mean-field level corresponds to a spin 
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nematic [ pl[ . Indeed, let us first introduce the z-component of tlie spin current J^^ associated 
to tlie ath spin-S XY cliain (a — 1, 2): 

JL^-v\l^d,Qa+. (31) 

Tfie vacuum expectation value of this operator can then be computed since one has in the 
ground state of the mean-field Hamiltonian (|27|): {dxO+) = —k/v ^ and (dxO-) — 0. 
This latter result stems from the fact that the Hamiltonian (7i-) in Eq. ( p7| ) is a standard 
sine Gordon model characterized by a ground state with zero topological charge. Using the 
redefinition (p6|), one finally obtains the following estimate: 

(JL) - {JD - {d.e+) = ^ 0. (32) 

These spin currents can also be expressed in terms of the original spin degrees of freedom of 
the lattice Hamiltonian ( p3| ) using the identification (p^): 

((4nA4„+i)j ~-A2y^(a,e+) ^0, a = 1,2, (33) 

whereas similarily the (interchain) zigzag spin current along the z-axis reads as follows: 

' S , . ( [tt - \, ^ S 



Ji{[Si,n A S2.n)j ^ -2^-g(sin (^^ I©-]) = ^ (34) 

where Eq. ( ^ ) has been used. 

The physical picture that emerges from this mean-field analysis is therefore a spin nematic 
phase that preserves the U(l) and time reversal symmetries and displays long-range chiral 
ordering in its ground state (^, . In the classification of Ref . , this phase corresponds 
to a p-type spin nematic. At this point, it is important to stress that this chiral ordering 
is different from the scalar chirality order operator |32j: {Si^n- (^S2,n A S2.„-i^) which breaks 
parity and time reversal symmetries. In our case, the spin nematic phase does not break the 
time reversal symmetry but spontaneously breaks a Z2 symmetry of the model which, as it will 
be shown in the next section, is a tensor product of a site-parity and link-parity symmetries on 
the two chains. As a result, as first discovered in the S=l/2 case in Ref. |3j, this produces a 
picture of local nonzero spin currents (|3^, polarized along the z-anistropy axis circulating 
around the triangular plaquettes of the two-leg zigzag spin ladder. 



B. stability of the chiral critical phase 



It is important to study further the stability of this critical spin nematic phase (chiral critical 
phase) in the + channel with respect to various operators that will be generated in higher orders 
of perturbation theory or equivalently terms consistent with the symmetries of the original 
lattice model. Indeed, on general grounds, one expects that some operators in the -I- sector 
should destroy the criticality of the phase at least for some finite value of Ji/J2- First of all, 
in the mean-field approach, the twist term acts like a sort of magnetic field. As is well known, 
a magnetic field along the anisotropy axis is a source of incommensuration but also leads to 
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a renormalization of the compactification radius of the bosonic field. This last effect was not 
found in the previous approach as seen in the universal behavior of the spin-spin correlations 
(pO|). From a symmetry point of view (continuous U(l) symmetry), there are no reasons to 
expect such universal behavior. It is a first sign that higher order terms in perturbation theory 
might be important here. On the other hand, as seen in Section II, there is at least a massive 
region (the dimerized or Haldane phases) in the phase diagram when increasing the value of Ji 
at fixed J2- It is therefore likely that a vertex operator, generated in the renormalization group 
flow, in the + channel will kill the critical phase at least for a critical value {Ji/J2)c- 

We shall now discuss the bosonic representation of the different discrete lattice symmetries 
of the model ( p3| ) to find the nature of the operator that will be generated in the + sector by 
the renormalization group flow. Let us flrst consider the one-step translation (t^aj ) , site parity 
(Pg"^), and link parity (P^^) corresponding to the chain of index a — 1,2. Using the deflnition 
(0) of the diagonal bosonic fleld that accounts for the criticality of the XYl spin fluid phase 



in the decoupling limit (Ji = 0) and the bosonic representations ( [Aq , AE , AlC ) in the S=l/2 
case described in the Appendix A, one obtains the following identifications respectively for the 
one-step translation, site parity, and link parity: 



Qa+ e>a+ + + p^VS^, (35) 



ea+ [X) ea+ i-x) + q'aVS^, (36) 



and 



$a+ {x) -$a+ i-x) + 7la\j ^ 

Qa+ (x) ea+ i-x) + V2^ + n^Vs^, (37) 

where Pa,Pa,Qa,qa^''^a, are integers. From these correspondences, one can deduce the bosonic 
representations of the discrete symmetries of the Hamiltonian (p3|). The translation symmetry 
acts on the symmetric and antisymmetric combinations ( p6| ) of the bosonic fields <I>q+ as follows: 

4-+ ^ 4-+ + + (pi + P2) 

e+ ^ 6+ + V45^ + Va^ (p'l + P2 

i-^$- + y^(pi-P2) 
e_ ^ e_ + \/4^ (p'l - • (38) 



A second type of discrete symmetry of the Hamiltonian (|2^) si consists of a vertical axial 
symmetry combin 
in the following): 



symmetry combined by an one-step translation symmetry t'^aa along the lower chain (labelled 1 
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S2.n — * S; 



2,-n, 



(39) 



namely in the continuum limit: 



fii {x) — ni (— x) 
^2 (a;) n2 {-x) , 



(40) 



which corresponds to a tensor product of a hnk-parity transformation on chain 1 and a site- 
parity transformation on chain 2 (si = P^^ ® P^"^) when the model is viewed as a zigzag ladder 
(Eq. (po|)). The bosonic representation of this discrete symmetry is thus: 



(^) 

0+ i^) 

(x) 

e_ (x) 



'7r5 



(?^i + 92) 



' 2 V 45 
e+ (-.x) + \/5^ + Vi^ (n[ + (72) 

(~x) - ^- + y ^ ("1 - 92) 



e_(-x) + \/SV + \/4^(?i'i-92) ■ (41) 
In the same way, the Hamiltonian (^3) is also invariant under the transformation (s2 symmetry): 



51. n 

52, n 



Sl.-n 



(42) 



which can be viewed as a P^^-* ® P^"'^-' transformation. In terms of the bosonic fields of the basis 



, this latter symmetry is realized through: 



(x) 

e+ (x) 
!>_ (x) 
e_(x) 



^$+(-x) + 



tS 



{n2 + qi) 



e+ (-x) + ^fs^ + [n^ + 91) 



(-x) + 

e_(-x)-v 



(-712 + qi) 



2 ■ V 45 
57r + V 47r5 ^— ri2 + 



(43) 



There is a second family of discrete symmetries of the Hamiltonian (|2^): S3 — P12 <8i t'^o or 

S4 = P] 



12 ( 



(2) 



which corresponds to an interchange of the chains combined with a translation 
symmetry along the lower or upper chain. In terms of the original spin degrees of freedom, the 
S3 and S4 symmetries respectively write: 



Sl^n 52, n 
52, n —> 5l,n+l, 



(44) 



51, „ 

52, n 



52,n-l 
5l,n, 



(45) 



so that in the continuum limit, one has 
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and 



Hi (x) - 
"2 (a;) - 

ni (x) 
n2 (x) 



^2 (x) 
-ni {x) , 

-n2 (x) 
ni (x) . 



(46) 



(47) 



The bosonic representation of these last discrete symmetries of Eq. (|2^) is then respectively 
given by: 



e+ 



$+ + 

e+ + 



%5 



-Pi 



2 V 45' 



(48) 



Oh 



$+ + 
e+ + 



'ttS 



P2 



2 ■ V 45 
^ + \/47r5p2 



45^^ 



'SV + v47r5p2' 



(49) 



With all these identifications, one observes that the continuum limit (|2^) of the lattice Hamil- 
tonian (23) is invariant under all discrete symmetries (|3^, ^9|) as it should be. 

However, the mean- field Hamiltonian ( p7| ) is invariant under (|3^, and but breaks the 
si,2 symmetries ( ^H |4^ ). These latter Z2 discrete symmetries are spontaneously broken in 
the ground state of the critical spin nematic phase and account for the formation of nonzero 
local spin currents polarized along the z-axis circulating around the triangular plaquette of 
the two-leg zigzag spin ladder. The operator, that occurs in the + sector of the mean-field 
Hamiltonian (|27|), with the smallest scaling dimension and consistent with the symmetries (38, 
|4^ , ^ ) without breaking the continuous U(l) diagonal symmetry of the model turns out to be: 
cos(/i\/167r5i'+), with ^ = 1 (respectively ^ = 2) if S is integer (respectively half-integer). The 
stable effective field theory in the + channel is thus: 



Kd:,Q+ - ^ cos 



5$. 



(50) 



where the value of the Luttinger parameters f , K cannot be determined within this bosonization 
approach. For a small value of J1/J2 (i.e. K ~ 1), the cosine operator in Eq. ( |5C| ) is a 
strongly irrelevant contribution and the system displays a critical phase with incommensuration 
generated by the dxQ+ field. This chiral critical phase, first predicted in the S=l/2 case in Ref. 
1^, is thus a generic phase in the large J1/J2 limit of the model (|^) in the general spin-S case. 
In particular, it is worth stressing that, in the S=l/2 case, the operator cos (•\/87r$+), which 
opens a mass gap in the + channel and thus destroys the chiral critical phase found in Ref. 
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is not generated by the renormalization group flow. Indeed, while this latter operator is 
permitted by the translation symmetry (^8|), it is odd under the S3 and S4 discrete symmetries 



( [4q , 49) which forbid its presence in the low-energy effective field theory. This result leads 
us to expect that the chiral critical phase does exist in the certain range of the parameter of 
the lattice model for S=l/2 in full agreement with the very recent numerical study As 
J1/J2 is further increased, it is natural to expect that the effective theory ( |50| ) describes a 
phase transition of KT type from the chiral gapless phase at geff = to a chiral gapped phase. 
Indeed, there will be a critical value (Ji/J2)c (the Luttinger parameter at the transition being 
equal to Kc — l/(2S'/i^)), which cannot be obtained within this bosonization approach, above 
which the cosine operator cos(/3$+) becomes relevant and a mass gap opens in the + sector 
(KT transition) without killing the incommensuration stemming from the dxQ+ operator. In 
this respect, this mechanism of generation of incommensuration is different from the usual 
commensurate-incommensurate scenario | ^3| since, in this latter case, there is a competition 
between the uniform spin density dx^+ field and the cosine operator cos(/3$+) leading to a 
threshold above which the incommensuration settles in the system. One should note that the 



existence of this incommensurate gapful phase when the cosine operator in Eq. (50) becomes 
relevant can also been seen using a Luther-Emery or Toulouse limit of the Hamiltonian (|5^) as 
it has been used to explain the origin of the incommensuration found in the phase diagram of 



the quantum axial next-nearest-neighbor Ising chain 1 34 1 . 

The full characterization of the intermediate phase (chiral gapped phase) depends on whether 
S is integer or half-integer. Indeed, for J1/J2 > (^1/^2)0, the bosonic field of Eq. ( pO| ) is 
locked in one of the minima of the potential —geff cos(^Vl67rS'l>+) which for geff > |Q are 
located at: ($+) = p ^ tt/AS/^, p being an integer. Moreover, the value of the compactification 
radius of the bosonic field is equal to: Rs — l/vl67r5. This follows from the redefinition 
( p^ and the fact that the compactification radius of the bosonic field that accounts for the 
critical properties of the spin-S XY chain is Rs — I/VSttS as it has been found in Section II. 
From the precise knowledge of Rs, one deduces the following identification: 

^ #+ + 27ri?s = + (51) 

From this equivalence and the position of the minima corresponding to the pinning of the 
bosonic field we thus conclude that in the integer spin case {ji — \) the ground state of 
the massive phase is non-degenerate whereas for half-integer spins (/i = 2) there is a two-fold 
degenerate ground state j3|]. Therefore, the chiral gapful phase corresponds to a massive phase 
with a coexistence of incommensuration and a Haldane phase (respectively dimerized phase) in 
the integer (respectively half-integer) spin case. From the identification of the massive phase 
found at large J1/J2 in the weak coupling analysis (see Section II), we then expect an Ising 
(Z2) transition between the chiral gapped phase and the Haldane or dimerized phases as J1/J2 
is further increased. At this Ising critical point, the total spin current {dxQ+) vanishes i.e. 
the disappearance of the incommensurate behavior and the systems enters a commensurate 
massive phase: Haldane or dimerized phases depending on the spin. At this point, one has 
to mention that the existence of this intermediate incommensurate massive phase, within our 
mean-field approach, relies on the decouping of the degrees of freedom in the two channels 
+ and — as in Eq. (^). We cannot rule out a different scenario that might occur in the 
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system nonperturbatively due to the effect of the interactions in the two sectors: a single phase 
transition between the chiral critical phase and the Haldane or dimerized phases. At this 
critical point, one has simultaneously the appearance of a mass gap in the spectrum as well as 
the cancelation of the spin current so that the chiral gapful phase shrinks to zero in this case. 

IV. CONCLUDING REMARKS 

In the present work, we have investigated the low-energy physics of the one-dimensional spin- 
S Ji — J2 XY model within the bosonization approach. Around the two limits (J2/J1 1, 
J1/J2 ^ 1) where a field theoretical analysis can be performed, we have described the nature 
of the different phases that occurs as well as the determination of the effective field theories of 
the resulting phase transitions. The critical XYl spin fluid phase at J2 = is generically stable 
upon switching on a nonzero value of the next-nearest-neighbor interaction except for the very 
special S=l case where the Haldane phase is immediately stabilized in full agreement with the 
DMRG study of Refs. |^,^. As the exchange interaction J2 is further varied, the model exhibits 
a KT phase transition described by a standard sine Gordon model between the XYl spin fluid 
phase and a fully massive dimerized or Haldane phases depending on the value of the spin. In the 
zigzag ladder limit (J1/J2 ^ 1), we have shown that, whatever the value of the spin, the chiral 
critical phase, first predicted in the S=l/2 case by Nersesyan et al. ||^, should exist in a certain 
range of the parameters of the model. This interesting spin nematic phase preserves the U(l) 
and time-reversal symmetries but spontaneously breaks a Z2 symmetry {Pi}^®Ps^"^) resulting on 
the formation of nonzero local spin currents in the ground state polarized along the anisotropy 
z-axis. Futhermore, the transverse spin-spin correlation functions are incommensurate with a 
wave vector qs — tt/oo ~ ( Ji/ J2)*'^^'^'^~^'' and decay algebraically with the distance with an 
exponent 1/(85') obtained within the mean-field approach used here. As the interchain J1/J2 is 
further increased, one expects the existence of a KT phase transition between the chiral critical 
spin nematic phase and an incommensurate gapful phase (chiral gapped phase). In particular, 
the effective field theory corresponding to this transition has been determined in this work. The 
nature of this chiral gapped phase corresponds to a coexistence of incommcnsuration stemming 
from the presence of nonzero spin currents in the ground state and a Haldane or dimerized 
phases depending on whether the spin S is an integer or half-integer. We then expect an Ising 
phase transition associated to the disappearance of the spin current between the chiral gapped 
phase and the standard Haldane or dimerized phases. The phase diagram found in this work 
is consistent with the predictions of the large-S study of Kolezhuk ||ll[] except for the special 
S=l case where the XYl spin fluid phase shrinks to zero. It will be very interesting if some 
extended DMRG studies can be performed in the S > 1 case to further shed light on the 
physical properties of the model as well as the possibility to extract the Luttinger parameters 
of the effective field theory (|5^). The different phase transitions reported in this work could 
also be investigated by means of a level spectroscopy analysis as in the one-dimensional spin-S 
XXZ Hcisenbcrg model [||. 
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Note added 

When this work was completed, we became aware of a very recent work by Hikihara et al. 
p9[ | who have investigated the 8=1/2,3/2,2 Ji — J2 XY chain using a DMRG analysis. They 
have found that the chiral critical phase appears in a broad region of the phase diagram in 
the general spin-S case in agreement with our work. Futhermore, the prediction on the decay 
of the spin-spin correlation (1/(85')) in the chiral critical phase found within the bosonization 
approach has been numerically verified. Finally, for integer spins (S=l,2), the authors of Ref. 

have reported the existence of the chiral gapped phase in a very narrow region of the 
phase diagram whereas in the half-integer case (8=1/2,3/2) it has not been identified within 
the numerical precision of the work 139]. 



APPENDIX: A THE XY CHAIN IN THE CONTINUUM LIMIT 



In this Appendix, we shall recall some well known facts on the continuum limit of the XY 
chain to fix the notations that will be used throughout this paper. 

The Hamiltonian of the antiferromagnetic spin-1/2 XY chain is (Ji > 0): 

Wo = Jl E + ^nS'n+l) , (Al) 

n 

where is a spin-1/2 operator at site n. As is well known, this model can be written in terms 
of lattice fermions Cn using the Jordan- Wigner transformation: 




(A2) 



The continuum limit of the model (Al) can then be performed with the introduction of right 
and left-moving fermion fields R,L: Cn/y/oQ — > R{x){iY^°'° + L{x){—i)^/°-°,x = nao, ao being 
the lattice spacing. Using the fermion-boson correspondence (see for instance Refs. js^,^): 

R = ^ exp (iV^^it) 
v27rao \ ^ 

L = exp (-i%/4^$L) , (A3) 

the Hamiltonian ( |Al[ ) can be expressed in terms of a bosonic field $ and its dual field O in the 
continuum limit: 

no = ^Jdx (^{d.ef + (9,$)") , (A4) 

where vo = Jiao is the spin velocity and we work with the following conventions: 

4> = <I>L + 

e = $L - 

[*«,$l]=«/4. (A5) 
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This latter commutation relation is necessary to insure the anticommutation between the right 



and left fermion operators (see Eq. (A3)). The bosonic field is compactified with the radius 
R = I/VItt: <f> ^ $ + ^/tF whereas the dual field is compactified with the radius R — l/{2nR): 
Q ^ + The spin density operator in the continuum limit decomposes into uniform and 

alternating parts: 

5~ J+(-l)"/"°n, (A6) 

which can also be expressed in terms of the bosonic fields as follows: 

-1 



■ sin I V 47r$ 

Trao 



exp (i^/wQ) 



= J- — (exp (i3V7r<I>L + i\/^$fi) + exp {~i3Vn<i>R - i\AF$L)) 



1 

: exp 



i^/¥e) sin (yi^^^ . (A7) 



-\/27rao 

We end this appendix by giving the bosonic representation of the discrete symmetries of the 



XY Hamiltonian (Al) that will be very useful when investigating the stability of the chiral crit- 
ical phase in Section III B. Under a one-step translation symmetry, the bosonic fields transform 
according to: 

Q^e + y/^ + pV^, (A8) 

p,p being integers since from Eq. ( |A6| ) the alternating part (n) of the spin density should 
be odd under the one-step translation symmetry. Under the site parity Ps {Sn S-n), the 
uniform and staggered parts of the spin density should be even so that: 

/tt 

$ {x) {-x) + ^+ iVt^ 



e{x) ^ e i-x) + q V47r, (A9) 

where q, q are integers. The link parity P^ (n — > 1 — n) is a combination of a site parity and a 
translation symmetry so that under P^ the bosonic fields $ and O transform as: 

$ {x) -$ {-x) + ny/n 

e{x) ^ e (-x) + y/n + n V4^, (AlO) 

n, n being integers. 



APPENDIX: B THE ONE-DIMENSIONAL S=l/2 Ji - Ja XY MODEL 

In this appendix, we derive the continuum limit of the model (|l|) in the S=l/2 case in the weak 
coupling limit J2 <C Ji- This calculation has been done several times j|,|l|,^ with different 
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bosonized expressions. This discrepancy stems from the fact that one has to be extremely 
careful when deriving the continuum limit and in particular for obtaining the correct velocity 
renormalization. We shall redo here this calculation for completeness and also since it will be 
needed in Section II when deriving the bosonization approach of the Ji — J2 spin-S XY chain 
in the J2 ^ Ji limit. 

The first step of the computation is to express the interacting part of Hamiltonian (|l|) in 
terms of the lattice fermions using the Jordan- Wigner transformation (|A^): 

n^nt = -J2 J2 (^1+2 (4+lCn+l - 1/2) C„ + i/.C.) . (Bl) 



Using the continuum limit of the fermions and the bosonization correspondence (A3) described 
in Appendix A, one has: 



Hi 



2tt 



dx ((-^)"/"'' : e-'^*« : [x + 2ao) 
+ : e'^^*^ : {x + 2ao)) 



( -^d^^ i.^ + ao) + ^ : sin (^V4^$j : {x + ao) j 

([if'"" : e'V^*H : {x) + : e"*^*^ : (a;)) + H.c. (B2) 

To derive the continuum expression of this Hamiltonian, we need the following operator product 
expansions in a standard gaussian c — \ theory: 



47r [z — w) v47r 



47r (z — w) v47r 

: (l-i^(z-w)d^R-2TTiz-wfid^RY\e'^'^^ : {w,w) 

2 [z — w) \ ^ ' ) 

■ e*v^*^ : (z) : sin fV4^$) : {w,w) "—^ : e*%/T6^*i(-)e^V47*«(») : 



+ r : f 1 + iVI^ (z - w) 9$L - 27r (z - w)^ (d^Lf) e-*^*« : («;, w) (B3) 

2 (z — ui) V / 

with the convention w — vqt + ix and dx = i{d — d). Using these results and keeping only 
non-oscillatory contributions in Eq. ( p32| ), we finally obtain: 

H,nt ^ / dxcos (\/l6^$) - / dx {d^e f (B4) 

which is in perfect agreement with the earlier derivation made by Haldane (see the erratum) 
and is in contradiction with some recent ones in the litterature [l|j3^. 
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